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Here we present neutron diffraction results on the mineral azurite. We have found that the crystal 
structure of azurite can be described in the space group P2i which is the next lower symmetric group 
of P2\/c as found in earlier works. This small change in symmetry does not greatly influence the 
lattice parameters or atomic fractional coordinates which are presented here for single crystal diffrac- 
tion refinements. The ordered magnetic moment structure of this material has been determined and 
is comprised of two inequivalent magnetic moments on copper sites of magnitude 0.684(14) and 
0.264(16) fiB- This result is discussed in terms of the anisotropic exchange and Dzyaloshinskii- 
Moriya interactions. It is found that the system is likely governed by one-dimensional behaviour 
despite the long-range ordered ground state. We also highlight the significance of strain in this 
material which is strongly coupled to the magnetism. 



PACS numbers: 75.25.-j 75.50.Ee 61.05.fm 75.30.Et 

I. INTRODUCTION 

Low-dimensional, quantum spin systems have gar- 
nered much attention of late from both experimental- 
ists and theorists alike. One such material is azurite, 
CU3 (003)2 (011)2, a natural mineral which has been pro- 
posed as an experimental realization of the ID distorted 
diamond chain model. Interest in the magnetic proper- 
ties of this material began in the 1950s when an ordered 
antiferromagnetic (AFM) ground state was observed be- 
low 1.9 K. 1 ^ From these studies, a detailed description of 
the magnetic moment structure was never revealed how- 
ever it has been speculated that the low ordering tem- 
perature may result from frustration in the diamond ar- 
rangement of spin- 1/2 Cu 2+ ions^. Early susceptibility 
measurements indicated that the coupling of spins within 
the chain could be described with alternating dimers and 
monomers 4 (Fig. [1} . This was later confirmed by magne- 
tization measurements which revealed a distinct plateau 
at 1/3 the saturation magnetization indicating a polar- 
ization of the monomer spins within an applied magnetic 
fields 

More recently, interest in azurite has been renewed 
with an ongoing debate about the relative strength of 
the exchange interactions in this material as defined in 
Fig. []|r.2. Current consensus is that this system can- 
not be described in terms of a simple isotropic exchange 
chain Hamiltonian, but rather that inter-chain coupling 
and anisotropic exchange must also be taken into account 
when describing this material. A detailed understanding 




FIG. 1: (Colour online) Diamond chain model showing the 
relative exchange interactions. It is now widely believed that 
J2 is the strongest interaction in azurite. 



of the magnetic structure, as presented here, may provide 
further useful material in the construction of an effective 
Hamiltonian to describe this system. 

To fully describe the magnetism in azurite, a complete 
understanding of the structural properties is also neces- 
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sary. A definitive set of lattice parameters and atomic 
positions should provide an accurate input for exchange 
coupling calculations. The crystal structure of azurite 
was already determined in the late 1950s by Gattow and 
ZeemaniS.. It was found that azurite crystallises in the 
monoclinic space group P2\/c (No. 14). This space 
group was confirmed by two independent studies by Zi- 
gan el alM- and Belokoneva et alr^. 

More recently, evidence of magnetoelastic coupling in 
this material has been revealed^. Confirmed by pre- 
liminary neutron scattering investigations, the struc- 
tural strain in this material coincides with the onset of 
magnetic ordering 14 . It is therefore possible that the 
structure of this material is appreciably different be- 
low Xj^ compared with the structure published for room 
temperature^. Given the interest in predicting the mag- 
netic properties of azurite from structural considerations 
alone it seems prudent to obtain accurate structural data 
on azurite in its magnetically ordered state. Here we 
present the results of a neutron diffraction study of azu- 
rite at temperatures above and below T^, revealing fur- 
ther evidence of structural distortion in addition to the 
ordered magnetic ground state. 



II. EXPERIMENTAL METHODS 

The sample used in this study was cut from a large 
high-quality crystal of azurite which has been used in pre- 
vious studies^. The sample was roughly cubic in shape 
(dimensions ~3x4x3 mm) and mounted on a copper 
pin which offered good thermal conductivity. Diffrac- 
tion measurements were performed on both the time-of- 
flight (TOF) instrument SXD at ISIS in the UK and the 
four-circle diffractometer D10 at Institute Laue-Langevin 
(ILL) in France. Additional diffraction measurements 
were also conducted using the El triple axis spectrometer 
(TAS) at Helmholtz Zentrum Berlin (HZB), Germany. 

The Single-Crystal Diffractometer SXD combines the 
time-sorted Laue method with a large array of position 
sensitive detectors to allow access to large volumes of 
reciprocal space in one simultaneous measurement. Us- 
ing this instrument, the sample is illuminated by a white 
beam of neutrons with a wavelength range 0.2-lO.OA. A 
full data set of Bragg peak intensities was measured at 
both 1.5 and 5 K for a comparison of the structural pa- 
rameters above and below the Neel temperature. 

Additional single crystal diffraction experiments were 
carried out on the instrument D10. Using a dilution 3 He 
cryostat, it was possible to determine structural and mag- 
netic properties at both 200 mK and 5 K. This ensured 
that the magnetic structure of the sample was deter- 
mined well below the ordering temperature of 1.87 K. 
D10 uses a Cu-monochromator selecting the single neu- 
tron wavelength A = 1.26 A. A selection of half-integer 
reflections and reflections forbidden by the P2i/c space 
group were also measured to ensure an accurate evalu- 
ation of the higher order scattering intensity and struc- 



tural symmetry respectively. The presence of multiple 
scattering could be excluded by the use of psi scans. The 
integration of the Bragg reflections was performed us- 
ing two different methods, the SEED program 15 and the 
RACER program 16 , in order to optimise the quality of 
data. The refinements of the crystal structure were car- 
ried out with the program Xtal 3.4i£. Here the nuclear 
scattering lengths 6(H) = 3.7409 fm, 6(H) = 6.6484 fm, 
6(0) = 5.805 fm, and 6(Cu) = 7.718 fm were usedi£. 

For the investigation of the magnetic structure of azu- 
rite, we collected data sets at D10 using the longer neu- 
tron wavelength A = 2.36 A. We collected a set of nu- 
clear reflections at 200mK in order to determine the over- 
all scale factor from the crystal structure refinements. 
With the absorption and extinction corrected magnetic 
structure factors, we were able to obtain the magnetic 
moments of the Cu 2+ -ions in the magnetically ordered 
range. The moments of the Cu-atoms were refined with 
the program FullProf 19 . The magnetic form factors of 
the Cu 2+ -ion was taken from Ref. 20. 

Finally to investigate the effects of applied magnetic 
fields on the structure of azurite, single-crystal neutron 
diffraction was carried out using the El triple axis spec- 
trometer at HZB (neutron wavelength A = 2. 428 A ). 
The small azurite single crystal remained on the copper 
mount with a horizontal a* — b* scattering plane, such 
that magnetic fields up to 14 T were applied perpen- 
dicular to b* , i.e., perpendicular to the chain direction. 
In this geometry the 1/3 magnetization plateau, where 
the monomer spins are polarised with the applied field, 
ranges from 11 to 30 Tesla 5 . A temperature dependence 
of these applied field effects was also conducted for tem- 
peratures up to 5 K. With the aid of a 3 He cryostat insert, 
base temperatures of T ~ 0.5 K were attained. 



III. RESULTS AND ANALYSIS 

A. TOF diffractometer measurements 

The data from SXD have revealed that Bragg peaks of 
index-type (hOl) with I odd were present in the scattering 
profile despite being forbidden by the P2\jc space group. 
This observation clearly indicates the loss of the c-glide 
plane of P2\/c. Group-subgroup relations between the 
space groups showed that Cu3(C03)2(OH) 2 crystallizes 
either in the noncentrosymmetric monoclinic space group 
P2i (No. 4) or the centrosymmetric triclinic one PI (No. 
2). For PI there is no limiting condition on any (hkl), but 
for P2i the type (OkO) is restricted by k being even. Due 
to the fact that reflections (OfcO) with k being odd could 
not be observed in our experiments, it can be concluded 
that the correct space group is probably P2\. The same 
limiting conditions on these reflections were observed at 
both 1.5 K and 5 K. This confirms the structural ori- 
gin of these peaks while indicating that azurite retains 
this symmetry both above and below the Neel temper- 
ature. The observation of these peaks with the TOF 
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method is unambiguous - the measured intensity cannot 
be attributed to higher-order scattering contamination, 
as it might with other diffraction techniques. In contrast, 
the 'forbidden' reflections were actually also observed in 
the room temperature single crystal neutron diffraction 
study of Zigan et at— but were attributed to higher-order 
scattering contamination. 

The atomic fractional coordinates were refined from 
the SXD data containing approximately 1800 unique 
structural reflections at the two measured temperatures 
of 1.5 and 5 K using the analysis program JANA2006 21 . 
In the analysis of the structural Bragg peak intensi- 
ties, information from data sets collected at different 
orientations of the crystal are not merged resulting in 
a large data set and complex wavelength-dependent ex- 
tinction corrections. This, coupled with additional wave- 
length dependent corrections, typically gives poorer Rp- 
factors for TOF methods when compared to monochro- 
matic techniques. The i?p-factor is defined as Rf = 
^2 \ \F a \ — \F C \\/J2 \F \ where F a and F c are the observed 
and calculated structure factors respectively. Refine- 
ments using the P2\/c space group gave a reliability of 
Rf ~ 0.080 and compare favourably with the powder 
refinement as seen in Table HI The refinement of the 
SXD data using the P2j space group gave unphysical 
anisotropic thermal parameters and as such did not give 
accurate positional parameters. 

The powder data of Cu3(C03)2(OH) 2 were also re- 
refined using the space group P2±. A comparison be- 
tween the refined lattice parameters from the P1\jc 
space group as published earlier— and the P2\ space 
group showed that despite the change in symmetry group 
the lattice parameters were the same within error. For 
the P2\ refinement of the powder neutron diffraction 
data taken at 1.28 K, the lattice parameters were a = 
4.99995(11) A b = 5.82256(14) A c = 10.33723(19) A 
with = 92.2103(17)°. 

Using the TOF technique with SXD, it was also pos- 
sible to observe and index the magnetic reflections be- 
low 1.8 K thereby confirming the propagation vector of 
k = (555), as has been observed in neutron diffraction 
measurements previously^. For absolute accuracy in the 
structure factors and somewhat more precise atomic pa- 
rameters further diffraction measurements at a continu- 
ous source were also conducted. 



B. Crystal structure from the four-circle 
diffraction measurements 

The crystal structure of Cu3(C03)2(OH)2 was refined 
using the data from D10 collected at both 200 mK and 
5 K. Refinements were carried out initially in the mono- 
clinic space group P2\/c as done earlier by Zigan et alr^ 
and Belokoneva et al*£. The lattice parameters were 
taken from the powder data refinement since they were 
considered to be more accurate. 
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FIG. 2: The temperature dependence of the integrated in- 
tensity of the (120), (040), (008) and (300) Bragg peaks of 
Cu3(C03)2(OH)2, as obtained from single- crystal neutron 
diffraction experiments. Note that the large increase in in- 
tensity at the Neel temperature is only observed for peaks 
with a non-zero Mndex, while peaks with a non-zero /i-index 
also exhibit a change in intensity below 0.5 K. Lines are a 
guide to the eye. 



For the refinement of the crystal structure at 200 mK a 
total number of 1771 (845 unique) reflections were used. 
The refinement of the overall scale factor, extinction pa- 
rameter and the positional and isotropic thermal param- 
eters resulted in a poor residual Rf = 0.089 (wRf = 
0.110). For the extinction correction in the structural 
refinement the formalism of Zachariasen (type I) was 
used 23 . The refineable parameter g in this formalism is 
related to the mosaic distribution and assumes a Gaus- 
sian distribution of mosaic blocks within the sample. An 
absorption coefficient /i — 0.158 mm" 1 , was also used 
for refinements. The refined extinction parameter (g) 
of about 1200 rad -1 indicates that extinction is quite 
strong in azurite. It could be seen that strong reflec- 
tions of the series (hOl) were observed to be systemati- 
cally weaker than the calculated values while the inten- 
sities of the (OfcO) reflections were calculated to be much 
stronger than their observations. This clearly highlights 
the anisotropy of the extinction. 

The significance of extinction effects in this material 
was revealed on measuring the temperature dependence 
of the Bragg peak intensities. As discussed in previous 
workis, such a large intensity increase below can- 
not be due to magnetic ordering alone. However since 
the intensity change coincides with the magnetic order- 
ing transition it is clear that the structure and magnetism 
of azurite are closely coupled. In Fig. f2J for example, we 
see a 30% increase in the intensity of the (120) peak on 
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TABLE I: Positional and thermal parameters of Cu3(CC>3)2(OH)2 as obtained from the structure refinements of single-crystal 
(on both D10 and SXD) and powder data in the monoclinic space groups P2\/c. The isotropic thermal parameters Ui S are 
given in units of 100 A 2 . For the refinement several thermal parameters were constrained to be equal. In these cases the 
standard deviation is listed only for one of the equal parameters. 
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X 
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U u 
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U u 


Cul 


2a 











0.40(4) 











0.14(1) 











0.12(3) 


Cu2 


4e 


0.2508(6) 


0.4967(5) 


0.0834(3) 


0.40 


0.2515(3) 


0.4976(1) 


0.08324(9) 


0.25(1) 


0.2516(2) 


0.4977(2) 


0.0834(1) 


0.12 


C 


4e 


0.3308(6) 


0.2994(5) 


0.3192(3) 


0.63(7) 


0.3294(3) 


0.2993(2) 


0.3180(1) 


0.26(2) 


0.3303(2) 


0.2992(2) 


0.3178(1) 


0.24(3) 


01 


4e 


0.0975(7) 


0.3972(6) 


0.3318(3) 


0.60(3) 


0.1012(4) 


0.3991(2) 


0.3310(1) 


0.36(2) 


0.1016(3) 


0.3992(2) 


0.3309(2) 


0.31(3) 


02 


4e 


0.0762(7) 


0.8126(6) 


0.4451(3) 


0.60 


0.0741(4) 


0.8126(2) 


0.4450(1) 


0.25(2) 


0.0743(3) 


0.8119(2) 


0.4450(2) 


0.28(4) 


03 


4e 


0.4518(6) 


0.2098(6) 


0.4183(3) 


0.60 


0.4500(4) 


0.2085(2) 


0.4176(1) 


0.35(2) 


0.4503(3) 


0.2093(2) 


0.4179(2) 


0.30(3) 


04 


4e 


0.4339(6) 


0.2949(6) 


0.2065(3) 


0.60 


0.4310(4) 


0.2937(1) 


0.2068(1) 


0.37(2) 


0.4317(3) 


0.2946(2) 


0.2072(2) 


0.37(4) 


H 


4e 


0.182(1) 


0.800(1) 


0.3709(7) 


1.75(15) 


0.187(1) 


0.7996(5) 


0.3682(3) 


1.72(6) 


0.1833(7) 


0.7996(6) 


0.3686(4) 


1.75(7) 



cooling the single crystal sample below Xj^ and an 80% 
increase in the intensity of the (040). On the other hand 
no change of intensity was observed for the (008). In 
fact only structural Bragg reflections (hkl) with a large k 
component were effected by a strong change of intensity. 
Surprisingly, additional yet weak intensity changes were 
also observed below 0.5 K on the (hOO) peaks suggesting 
that there may also be strain coupled to the a-direction. 
This appears at the same temperature as an anomaly ob- 
served in recent ultrasonic measurements^ however the 
nature of this anomaly is not yet known and requires 
further investigations. 

In order to improve the refinements we rejected 123 
(58 unique) of the very strong reflections. In Table Q] it 
can be seen that the positional and thermal parameters 
of the different atoms could be determined with good ac- 
curacy. The refinement finally resulted in a better resid- 
ual Rf = 0.070 {wRp — 0.069) although the positional 
parameters obtained from both refinements showed a 
good agreement. Further, the positional parameters of 
the single-crystal data also show a good agreement with 
the values obtained from the neutron powder diffraction 
data collected at 1.28 K and the single- crystal data from 
SXD (Table HJ. But despite the fact that the refinement 
quality from the single-crystal data was reduced due to 
anisotropic extinction we were able to determine the po- 
sitional parameters with much better accuracy than from 
the powder diffraction study (Table J}. The refinements 
of the D10 single-crystal data collected at 5 K resulted in 
a strongly enlarged extinction parameter of about 4000 
rad -1 which is clearly also affecting the fit parameters of 
the SXD data. 

It was again found in the DIO-experiment that reflec- 
tions of the type (hOl) with I being odd had weak yet 
significant intensity. The intensity profiles of the reflec- 
tion (101) at different temperatures is presented in Fig. 
[3J Finally a total of 60 unique reflections of this type 
could be collected. The presence of multiple scattering 
was checked by the use of psi-scans. Since the reflections 
of the type (OfcO) with k — odd were again not observed, 
it is likely that the crystal structure of azurite can be 
described in the space group P2\. However we carried 




Omega 

FIG. 3: Intensity profile of the [101] Bragg reflection at tem- 
peratures above and below Tj^. This reflection is forbidden 
for the P2\/c space group but allowed for P2\. Lines are 
guides to the eyes. 

out the structure refinements using both the monoclinic 
and triclinic space groups P2i and PI, respectively. The 
refinement of the single crystal data gave a residual Rf 
= 0.067 (wRf = 0.063). A significantly larger residual 
Rf = 0.076 (wRf = 0.074) was obtained from the refine- 
ment of 51 parameters using the monoclinic space group 
PI. A comparison of observed and calculated structure 
factors (in F 2 ) is shown in Fig. It was interesting 
to see that the calculated intensities of the reflections 
(010), (030), (050), (070) [(QA;0) with k = odd, forbid- 
den in P2i)] were practically equal to zero. This also 
indicates that the structure of azurite still contains the 
screw axis 2i parallel to b. Therefore we can conclude 
that azurite crystallizes in the monoclinic space groups 
P2i. The results of the refinement of the D10 data for 
the P2i symmetry group are summarised in Table [TT1 

The positional parameters of the copper atoms for the 
two different space groups P1\jc and P2i are given in 
Tables U and UU In P2\/c the Cul-atoms are located at 
the special Wyckoff position 2a(0, 0, 0). Therefore in this 
space group the interatomic distances to the next neigh- 
boring Cul '-atoms in (0,5,5), (0, — 5,5)) (0>5)~ f)> 
and (0,-5,-5) are identical [dcui-Cui 1 — 5.9321(1) A ]. 
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FIG. 4: Comparison of observed and calculated structure fac- 
tors (in F 2 ) for the P2\ (upper) and PI (lower) space groups. 
The distribution of points around the line F 2 bs — F 2 alc in- 
dicates the high quality of the fits. These data have been 
corrected for extinction effects as mentioned in the text. 



In the space group P2i the Cul-atoms are located at the 
general Wyckoff position 2a(x,y, z). In Table ITT1 it can 
be seen that the ^-parameter does not show any change 
from the ideal position x = 0. On the other hand the y- 
and z- values show a slight shift from the ideal position 
y = and z — j (in P2i) or z — (in P2\/c) of about 
3 — 4<7. In the lower symmetric space group P2\ one 
finds two different interatomic distances. The distances 
between Cul in (0, 0, 0) and the Cul'-atoms in (0, h, |), 
and (0, — |, i) are the same [e?i = 5.913(4) A ]. But along 
the opposite z-direction the distance to the Cul'-atoms 
in (0, and (0,-A,-|) is d 2 = 5.952(4) A where 

d\ and d 2 are plotted in the fee-plane in Fig. [5] In 
fact our study shows that the shift from the ideal po- 
sition is relatively weak. However it is interesting to note 
that the lower symmetric setting allows a dimerization of 
monomer Cul-sites. On the other hand the Cu2-atoms 
(in P2i/c) are located in the general Wyckoff position 
4e(x, y, z), while in the lower symmetric space group P2i 
this position splits into two different positions Cu21 and 
Cu22, but both at the general position 2a(x, y, z). In Ta- 
bles U and [XT] it can be seen that the Cu2-atoms in the 
two representations show a good agreement, with respect 
to their standard deviations. Please note that one finds 
the relation z'(Cu2) = z(Cu21 or Cu22) - \ for the two 
different representations. 



TABLE II: Positional and thermal parameters of 
Cu3(C03)2(OH)2 as obtained from the structure refine- 
ments of single- crystal data taken at 200 mK using the 
monoclinic space group P2i. For the refinement several 
thermal parameters, Ui S (given in 100 A 2 ), were constrained 
to be equal. In these cases the standard deviation is listed 
only for one of the equal parameters. 



Atom 


site 


X 


y 


z 


u is 


Cul 


2a 


0.0000(5) 


0.0015(6) 


0.2511(3) 


0.13(4) 


Cu21 


2a 


0.2514(4) 


0.4985(5) 


0.3336(2) 


0.16(3) 


Cu22 


2a 


0.7488(4) 


0.5027(5) 


0.1668(3) 


0.16 


Cll 


2a 


0.3306(5) 


0.3012(5) 


0.5685(3) 


0.21(3) 


C12 


2a 


0.6691(5) 


0.7026(5) 


0.9334(3) 


0.21 


Oil 


2a 


0.0997(5) 


0.3921(5) 


0.5815(4) 


0.27(3) 


012 


2a 


0.8954(5) 


0.5949(5) 


0.9198(4) 


0.27 


021 


2a 


0.0762(5) 


0.8100(5) 


0.6957(4) 


0.33(4) 


022 


2a 


0.9269(6) 


0.1853(5) 


0.8065(4) 


0.33 


031 


2a 


0.4518(5) 


0.2068(6) 


0.6700(4) 


0.33(3) 


032 


2a 


0.5518(5) 


0.7891(6) 


0.8341(4) 


0.33 


041 


2a 


0.4366(5) 


0.2915(5) 


0.4564(4) 


0.34(3) 


042 


2a 


0.5732(5) 


0.7037(5) 


0.0411(4) 


0.34 


Hll 


2a 


0.1768(12) 


0.8088(11) 


0.6129(7) 


1.60(7) 


H12 


2a 


0.8117(12) 


0.2113(23) 


0.8754(7) 


1.60(7) 



C. High magnetic field measurements 

Further insight into the magnetic and structural in- 
terdependence can be taken from the field dependence of 
the intensity of the (120) structural Bragg peak. For this, 
measurements were conducted using the El spectrometer 
at HZB. The left panel of Fig. [5] displays the integrated 
intensity of the (120) nuclear Bragg reflection as a func- 
tion of field while the right panel shows the temperature 
dependence. The intensity has been normalised either to 
the peak intensity at 14 T or at 2 K, since no further 
intensity changes were observed well within the plateau 
phase and above the ordering temperature respectively. 
At temperatures approaching T^, the field at which the 
intensity becomes invariant is reduced. 

The atomic fractional coordinates were not observed to 
change significantly within the accuracy of these neutron 
scattering experiments. However field and temperature 
induced changes in the lattice parameters have been ob- 
served in more accurate studiesi^^. These slight changes 
may facilitate the formation of the mosaic blocks within 
the crystal which in turn become misaligned on a macro- 
scopic scale increasing the mosaicity of the sample. 



D. Magnetic structure 

The magnetic order of the Cu 2+ -ions sets in at the 
Neel temperature Tj\j = 1.87 K. Previous diffraction 
experiments showed that the magnetic cell of azurite 
could be described with the propagation vector k — 
(1 11)^' ^ n or d er to determine the magnetic struc- 
ture of the Cu-sites we performed a symmetry analy- 
sis for the space groups P2\jc and P2i. This allowed 
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FIG. 5: (left) The field (left) and temperature (right) de- 
pendence of the normalized integrated intensity of the (120) 
Bragg peak of Cu3(CC>3)2(OH)2 near the plateau field range 
from single crystal neutron diffraction. Lines are guides to 
the eye. 



us to deduce the possible spin configurations compati- 
ble with the symmetry of the crystal structure accord- 
ing to the irreducible representations (ireps) T v as de- 
scribed elsewhere 2 ^. In a magnetic structure the mag- 
netic moment /i of an atom j at the lattice point R„ 
in the unit cell is characterized by the propagation vec- 
tor k, where the moment is given by the Fourier series 



-27rik-R,n 



The Fourier series Sj^ are 



linear combinations of the basis functions of the ireps 
obtained from the symmetry analysis. In the mono- 
clinic space group P2\/c the two different copper atoms 
are located (within the unit cell) at the Wyckoff posi- 
tion 2a[Cull(0, 0, 0), Cul2(0, |, I)] and 4e[Cu21(a;, y, z), 
Cu22(l -x,y + h, z), Cu23(l - x, 1 - y, 1 - z), and 



Cu24(x, -y + \,z+ 1)], with x = 0.2513(3), y = 

0. 4976(2), z = 0.0834(2). There are four symmetry oper- 
ators in P2\/c leading to four one-dimensional complex 
irreducible representations namely T\, r 2 , T 3 and T4. 
Since the propagation vector is k = (555) then k and 
-k are equivalent such that the Fourier series Sj^ should 
be real. The induced representation of the 2a site can be 
decomposed as r 2a = 3Ti + 3r 3 . The Fourier coefficients 
for the 2a site are then a linear combination of the irre- 
ducible representation Fi and its complex conjugate T 3 

1. e. (1 + i)Ti + (1 — i)T 3 leading to the following mag- 
netic configuration with real Fourier components, = 
(u,v,w) for the atom Cull and Sj tl2 = [u, — v, w) for 
the atom Cul2. The induced representation of the 4e 
site can be written as = 3Ti + 3T 2 + 3T 3 + 2>T^. 
In a similar way we obtain the Fourier coefficients for 
the 4e. In the lower symmetric space group P2\, Cul 
does not change its multiplicity and it is located at the 
position 2a[Cull(a;, y, z) and Cul2(l — x, y + 5, 1 — z)] 
with x = 0.0000(3), y = 0.0015(6), z = 0.2511(3). For 
the atoms Cull and Cul2 the Fourier components are 
found to be the same for both space groups P2\ and 
P2\jc. Due to the fact that the centre of symmetry 



is lost in P2\ the four Cu2 atoms split in two differ- 
ent sites, both in the same Wyckoff position 2a as also 
found for Cul. Within the unit cell the Cu2 atom are lo- 
cated in: Cu21(x, y, z), Cu22(l — x, y + ^, 1 — z) with x 
= 0.2514(4), y = 0.4985(5), z = 0.3336(2); Cu23(x, y,z) 
and Cu24(l - x, y - |, 1 - z) with x = 0.7488(4), y = 
0.5027(5), z = 0.1668(3). The only difference between the 
different setting is that the components u, v, w of the sub- 
sets (Cu21,Cu22) and (Cu23,Cu24) must not be necessar- 
ily the same. However, our crystal structure refinements 
showed that the positional parameters are only slightly 
changed in the lower symmetric space group P2\. There- 
fore the values of u, v, w were constrained to be equal for 
all the four Cu2-atoms. Further it is interesting to see 
that coefficients for the atoms Cul2, Cu22, and Cu24 are 
purely imaginary. Then a phase factor of j(in fractions 
of 2-71") between the atoms of one subset has to be taken 
into account. 

For the refinements of the magnetic structure we used 
a total number of 40 independent magnetic reflections. 
It could be shown that the coupling between the Cu- 
atoms of the monomers (Cull and Cul2) is compatible 
with spin arrangement of irep Ti: = (u,v,w) and 

Sj {12 = (u,—v,w). The moments are coupled parallel 
in the monoclinic ac-plane, and antiparallel along the b- 
axis. For subsets of the Cu2-site we found the spin con- 
figurations: Sj c21 = (u, v, w), = (u,v,w),S^ 23 



( 



-v, —w), S 



k24 



(- 



The spin coupling 



between Cu21 and Cu22, as well as the spin coupling 
between Cu23 and Cu24 along the monoclinic 6-axis is 
compatible with the symmetry analysis. 

The refined magnetic structure, resulted in a residual 
Rf = 0.069. Fig. [6] shows the refined magnetic struc- 
ture where the blue (dark) spins represent the Cu(l) sites 
which we commonly refer to as the monomer sites and 
the red (light) spins represent the Cu(2) sites which form 
dimerised pairs. Along the chain direction, the monomer 
spins are collinear and oriented in the ac-plane with a 
slight canting along the 6-direction. The monomer spins 
lie with an easy axis at an angle of 55 ±3° from the c-axes 
which corresponds within error to the value of 52° found 
from magnetisation measurements^. 

The magnetic moment on the monomer sites, is 
0.684(14) /is, while the moment on the dimer sites is 
0.264(16) fis- The relative x, y and z components of each 
site are shown in Table [TTTl A moment size of 0.684 fig is 
consistent with magnetic moments found in other cuprate 
systems where the reduction from the 1 [Ib expected for 
a Cu 2+ ion can be attributed to zero point spin fluctu- 
ations and/or covalency effect s 25 ' 26 . Note that the size 
of the moment on the dimers determined in this study, 
while significantly reduced from the monomer value, is 
somewhat larger than the moment found for azurite in 
its plateau phase from NMR-. This reduction may be 
caused by mean field effects. In the plateau phase, the 
monomers form a ferromagnetic chain which may influ- 
ence the overall moment at the dimer sites. In low fields, 
when the monomers are coupled antiferromagnetically, 



7 



TABLE III: Table of the magnetic spin components of the 
monomer (Cu(l)) and dimer (Cu(2)) spins at the atomic frac- 
tional coordinate positions of Cu(l) and Cu(2). 



Site 


Spin-component (/is) 


Full moment ((j,b) 


Cu(l)« 


0.42 (2) 




Cu(l) v 


0.36 (3) 




Cu(l), 


0.42 (3) 


0.684 (14) 


Gu(2) x 


0.036 (17) 




Cu(2)„ 


0.243 (15) 




Cu(2) z 


0.099 (16) 


0.264 (16) 



the mean field effects are somewhat different. 




b 



FIG. 6: (Colour online) The magnetic structure of azurite, 
shown in the b-c, a-c and a-b planes (l.-r.). The monomer 
sites, with a moment size of 0.684 /is, are represented by the 
blue (dark) spheres, while the dimer sites (0.264 ^b) are rep- 
resented by the red (light). The monomer spins on alternate 
a — b planes are aligned perpendicularly. The lengths d\ and 
d,2 have been determined for the P2\ space group in section 
C. 



IV. DISCUSSION 



Extinction 



Due to the close coupling between the extinction and 
the onset of magnetic order, it is likely that the magne- 
toelastic coupling is mediated by the diamond chain units 
which propagate along the crystallographic &-axis parallel 
to (OfcO). As a magnetic field is applied perpendicular to 
the 6-axis, the AFM correlations between monomer spins 
gives way to ferromagnetic correlations which in turn al- 
lows the lattice to return to a state with fewer mosaic 
block discontinuities and therefore increased extinction 
effects. At temperatures closer to T^, the magnetic in- 
teractions are weaker and thus a lower applied field is 
required to reduce the spread of mosaic blocks. This 
may explain why we see a reduction in intensity of the 
nuclear (120) Bragg peak in applied fields as seen in Fig. 
[7] Thus the change in intensity with applied field, whilst 
being an extinction effect, highlights the phase boundary 
of the AFM state in azurite. It is also worthwhile to note 
that the intensity of the (120) Bragg peak in the plateau 
phase is the same as in the paramagnetic, low field phase 
indicating that the mosaicity in both phases is the same. 



15 



10 
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1.0 1.5 

T(K) 

FIG. 7: An indirect measure of the phase diagram as inferred 
from the field and temperature dependence of the integrated 
intensity of the (120) peak (from Fig. [5}. Line is a guide to 
the eye. Ovals and arrows represent the Cu-dimer pairs and 
the relative orientation of the Cu-monomers repsectively. 



B. Magnetism 



All magnetic Bragg peaks showed an increase in in- 
tensity below the ordering temperature as expected for 
a second order phase transition. Since the influence of 
extinction is roughly proportional to the intensity of the 
Bragg peaks, the effect of the anisotropic extinction on 
the overall magnetic refinement result is minimal. 



For a reduced moment size of 0.264 /is we can con- 
sider that the Cu(2) sites (represented by the red spins 
in Fig. [SJ are coupled primarily into spin singlet states 
with their nearest neighbouring Cu(2) spin in the dia- 
mond unit. The observation of a finite spin polarisation 
implies that the singlet state is subject to some perturba- 
tive effect which induces the non- negligible spin moment. 
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Assuming the interchain coupling is negligible in compar- 
ison with the intrachain coupling, the staggered field of 
the neighbouring monomer spin sites may be responsible 
for this perturbation. 

Considering at first only Heisenberg exchange cou- 
plings, the expectation values of the spin moment on the 
dimer atoms in the staggered field of the neighbouring 
monomer sites can be calculated from the Hamiltonian 

H = J2S1.S2 + ( J\ — J3) (S m ) S* ~ (yi — ^3) ("S'm) $2 (1) 

where the ordered moment on the monomer sites de- 
fines the z-axis and is of magnitude (S). For all nonzero 
(S) a staggered moment on sites Si and £2 is established 
for the ground state collinear with the monomer spins. 
Whilst the ratio of (Ji - J 3 ) / J 2 = 1/2.8 implied by the 
relative ordered moments on dimer and monomer spins 
is in rather good agreement with the exchange couplings 
determined by Rule et al£ the non-collinearity requires 
the introduction of other exchange terms. 

A likely scenario involves Dzyaloshinskii-Moriya (DM) 
interactions which strongly influence the non cen- 
trosymmetric J\ and J3 exchange interactions as 
outlined in Figure [5] Taking the extra ex- 
change terms of the type D x . (§ ml X Si + S m 2 X S2 I ~t~ 

D3. (^S m 2 x Si + S m i x £2) and using the fact that 

S m i = —S m 2 = (0,0, (S m )) then the Hamiltonian be- 
comes 

H = J 2 S 1 .S 2 + (Df - D\) (S* m ) (Sf - 5f) 
- (D*i-D%)(S* m )(Sl-S%) 
+ (Ji - Js) (S z m ) (SI - Si) 

This then causes a titling of the dimer spins compared 
to the monomer z direction along the new effective stag- 
gered field (Df - Df , Df - D% , J x - J 3 ) (S z m ) direction. 
The significant tilting observed would imply that the DM 
interactions are of order (Ji — J3), i.e. a few K. 




FIG. 8: (Colour online) Dzyaloshinskii-Moriya (DM) interac- 
tions on the diamond shape units where the vertical direction 
corresponds to the fe-axis. We define the notation of and ® 
to represent antiparallel vectors since the actual orientation 
of the DM vector is not known. 

A further consequence of including DM interactions 
is that the effective Hamiltonian of the monomer spins 



loses full rotational symmetry and results in an XXZ type 
Hamiltonian with planar character (the easy plane be- 
ing perpendicular to (Di — D3)). To explain the spin 
flop type behaviour observed in azurite 3 ^ other terms re- 
sulting in an easy axis are needed, either from the sym- 
metry lowering to P2i, interchain coupling, or further 
anisotropies. The application of non-commuting mag- 
netic fields to such low dimensional XYZ antiferromag- 
nets, where J x ^ J y ^ J z , have been observed to cause 
remarkable behaviour in the vicinity of quantum critical 
point o 31 ' 32 and azurite may well provide an important 
new system in this line. Here, for example, the quantum 
phase transition into the plateau phase could take on a 
transverse Ising character—. We defer further discussion 
of the full derivation of the Hamiltonian to a future paper 
on the dynamics in azurite. 

The question remains as to what extent the diamond 
chains in azurite can be considered to be modeled reason- 
ably as a one-dimensional system. The Mermin- Wagner 
theorem indicates that truly one-dimensional systems 
cannot enter a long ranged ordered state due to quantum 
fluctuations 2 ^,, however experimentally, materials such as 
azurite which display quasi-lD characteristics clearly en- 
ter a 3D Neel state. The significance of quantum fluctua- 
tions in a given system may be established from the ratio 
of the Curie- Weiss temperature to the Neel temperature. 
It has been found from susceptibility measurements that 
for azurite this value is roughly —Qcw/Tn ~ 5, imply- 
ing that quantum fluctuations should play a significant 
role for all temperatures below 6cw^- 

Further, the relative strengths of inter- and intra-chain 
couplings may be inferred from the ordered magnetic mo- 
ment for comparatively weakly coupled chain a 3 ' 35 . The 
effective one-dimensional chain model provides a reason- 
able description of the system with coupling J mono = 
0.87 meV, determined previously from inelastic neutron 
scattering!. Taking the measured moment of too — 
0.684/is, then an effective tetragonal interchain coupling 
of |J_i_| = 0.097 meV would be required to account for 
the ordering strength. It should be noted that Ising 
anisotropies would serve to increase the ordered moment 
for the same interchain coupling. 



V. CONCLUSION 

In summary, we have observed that azurite belongs 
to the lower symmetry space group P2\. We have re- 
vealed the significance of magnetoelastic strain in this 
material as observed by field and temperature depen- 
dent extinction effects. Also presented here is the ground 
state magnetic structure of azurite which confirms that 
the diamond chain arrangement of magnetic Cu 2+ sites 
in this material can, to some extent, be considered as an 
alternating arrangement of dimer and monomer entities. 
The ordered magnetic structure may also indicate the 
presence of competing interactions between the chains in 
this material and the significance of anisotropic exchange. 
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Comparison of the magnitude of the ordered magnetic 
moment with existing theories of coupled quantum spin 
chains implies that the interchain coupling is weak in 
this material and that the system may thus be consid- 
ered one-dimensional. The magnitude and orientation 
of the magnetic moments in azurite points to additional 
anisotropy terms in the Hamiltonian. Additional diffrac- 
tion studies will be required to determine the ordered 
magnetic state of this intriguing material in its plateau 
phase. 
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